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Abstract 

The generalized Gauss-Bonnet theory, introduced by Lagrangian F(R, G), 
has been considered as a general modified gravity for explanation of the 
dark energy. G is the Gauss-Bonnet invariant. For this model, we seek 
the situations under which the late-time behavior of the theory is the de- 
Sitter space-time. This is done by studying the two-dimensional phase 
space of this theory, i.e. the R — H plane. By obtaining the conditions 
under which the de-Sitter space-time is the stable attractor of this theory, 
several aspects of this problem have been investigated. It has been shown 
that there exist at least two classes of stable attractors: the singularities 
of the F(R,G), and the cases in which the model has a critical curve, 
instead of critical points. This curve is R = 12H 2 in R-H plane. Several 
examples, including their numerical calculations, have been discussed. 

1 Introduction 

Based on various observations, it is believed that our universe is now in an 
accelerating phase [1]. Although the origin of this accelerated expansion is not 
yet known, almost all data indicate that nearly 70% of the present universe is 
composed of dark energy, the physical object that induces the negative pressure. 

There are two main classes of models that have been introduced as candidates 
of dark energy. The first class is based on the Einstein cosmology but with extra 
physical object as the source of dark energy. The scalar field (one-component 
or multicomponents) models [2], the scalar-tensor theories [3] and the k-essence 
models [4] are examples in this context. 

The second class of the models is based on the assumption that the gravity 
is being (nowadays) modified. The simplest one is obtained by adding a cos- 
mological constant term to Einstein action. This model suffers two important 
problems known as cosmological constant and coincidence problems [5]. Also, 
the cosmological constant model is a static model of dark energy and has not 
any dynamical behavior. The other modified gravity models are those that are 
based on the new actions. 
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The first family of these modified gravity theories are those known as f(R) 
gravity, with the action 



S = / d 4 x^^ 



(1) 



In h = c = G = 1 units, k 2 = 87r, R is the Ricci scalar and C m is the Lagrangian 
density of dust-like matter. Many features of f(R) gravity models, such as local 
gravity tests, have been studied [6]. 

Another well motivated curvature invariant, beyond the Ricci scalar, is the 
Gauss-Bonnet (GB) term 



G — R 2 — AR^vR^" + R^y^crR^^ ') 



(2) 



which is inspired by string theory [7] and is a topological invariant in four 
dimensions. The second family of modified gravity theories is known as f(G) 
gravity and is defined through 



S = 



d A x, 



-fJ 



^L R +f(G) + £„ 



(3) 



This model has gained special interest in cosmology [8] and its coupling to 
scalar fields, as it naturally appears in low-energy string effective actions [7], 
can introduce extra dynamics to this model. Other aspects of modified GB 
gravity, such as the possibility of describing the inflationary era, transition from 
deceleration phase to acceleration phase, crossing the phantom-dividc-line, and 
passing the solar system test have been discussed in [9]. 

The natural generalization of action ([3]) is the generalized Gauss-Bonnet 
dark energy, which have been introduced in [7, 10] : 



S = 



[F(R,G)+jC r , 



(4) 



Clearly the f(R) gravity and /(G) gravity are special examples of modified 
F(R, G) gravity. The hierarchy problem of particle physics and the late-time 
cosmology have been studied in F(R, G) models [11]. Recently, the behavior of 
these models in phantom divide line crossing and deceleration to acceleration 
transition, including the contribution of quantum effects to those phenomena, 
have been studied in [12]. It has been shown that the quantum effects can 
induce these transitions, when they are classically forbidden. 

One of the important characteristics of all dynamical systems, including 
the dynamical models of dark energy, is their late-time behaviors, studied in a 
framework known as attractor solution of dynamical systems, which has been 
deeply studied in mathematics. For dark energy models, the attractor solutions 
of scalar theories and some of the modified gravity theories have been studied 
in [13-15]. 

The main step in studying the attractor solution of a dynamical system is 
considering a set of suitable dynamical variables x%(t), . . . , x n (t), such that their 
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first time derivatives dxi/dt do not depend explicitly on time: 

— = Fi(xi, . . .,x n ), 

: (5) 

dx n . , 

n , " n\X\, • ■ • , X n ) . 

at 

The space constructed by variables xi,...,x n is called the phase space of the 
system and the system of equations ((5)) is said to be autonomous. The functions 
Xi (t), . . . , x n (t) define a path in the phase space and there is a unique path which 
passes any specific initial values xi(io) — %i,q> ■ ■ ■ ,x n (to) — x n .o, i.e. the paths 
do not intersect one another. The only exception to this statement occur at 
points (xx <c , . . . , x rhC ) where 

Fi(a;i !C , . . . , x ntC ) = 0, . . . , F n (x llC , . . . , x n>c ) = 0. (6) 

These points are called critical points, and any paths near these points, under 
specific conditions, will lead them at t — ► oo. In these cases, the critical points 
are called the stable attractors. 

The present paper is devoted to the study of the phase space and attractor 
solutions of generalized Gauss-Bonnet dark energy models. We will consider 
the R — H space (R is the Ricci scalar and H is the Hubble parameter) as the 
phase space of these models and show that in special case of F(R, G) = f(R), 
the results of [16], in which some features of phase space of f(R) have been 
studied, are reproduced. The choice of this phase space, which is the only 
possible choice in general F(R, G) model, has an important property. Since the 
attractor solutions are those which asymptotically lead to R = and H = 0, or 
R = R c and H = H c (R c and H c are some constant values), our phase-space 
study is in fact the study of possible de-Sitter solutions of modified generalized 
GB gravity. 

The scheme of the paper is as follows. In section 2, the set of autonomous 
equations of F(R, G) models is obtained and the condition of existence of stable 
attractors is discussed. Some specific examples of F(R, G) models which admit 
the stable attractors are investigated in section 3, and it is shown that the 
numerical studies confirm our results. Section 4 is devoted to the cases where the 
standard linear approximation method, used in obtaining the stability behavior 
of solutions, does not work. In section 5, it is shown that the singular points 
of the Lagrangian are always the stable attractors, and finally in section 6, the 
interesting cases where the critical points replaced by critical curves are studied. 
It is shown that all these critical curves always behave as stable attractor curves. 
We end the paper with a conclusion in section 7. 
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2 Critical points of F(R, G) gravity 

Consider the generalized GB dark energy model with action (j4]). Variation of 
this action with respect to the metric g^ v results in [11] 

~g>* v F(R, G) - 2F G (R, G)RR fiV + AF G (R, G)R p p R vp 

- 2F G (R, G)R ppaT R v prJT - 4F G (R, G)R wav R pa + 2(W p W l 'F G (R, G))R 

- 2g" u (W 2 F G (R, G))R - A{V P V P F G {R, G))R up - i(W p W u F G {R, G))R" P 

+ 4(V 2 F G (i?, G))BT + 4 3 ^(V p V CT F G (i?, G))R pa - 4(V p V CT F G (i?, G))R^ a 

- F R (R, G)R^ + WV U F R (R, G) - g pu W 2 F R (R, G) = 0. 

(7) 

Here, for simplicity, we do not consider the background matter field, i.e. C m — 0. 
In equation J?}, F R and F G are defined as following : 

. (8) 

For background metric, we consider, as usual, the spatially flat Friedmann- 
Robertson- Walker (FRW) metric in co-moving coordinates (t, x, y, z) as follows: 

d.s 2 = -dt 2 + a 2 (t){dx 2 + dy 2 + dz 2 ), (9) 

in which a(t) is the scale factor. The (t, i)-component of evolution equation ([7| 
then becomes 

-6H 2 F R (R, G) =F(R, G) - RF R (R, G) + 6HF R {R, G) 
+ 24H 3 F G (R,G)~GF G (R : G). 

H = a(t)/a(t) is the Hubble parameter. For this metric, the Ricci scalar R and 
the Gauss-Bonnet invariant G are 

R = 6(H + 2H 2 ), (11) 

and 

G = 2AH 2 {H + H 2 ), (12) 

respectively. The equations ((TX)|) - (fT2"]) are the Friedmann equations of F(R, G) 
gravity. The sum of (i, i) components of eq.(J7|) is obtained by using the time 
derivative of eq. (fT0|) and the eqs. lfTTj) and (fT2|) . 
From cq. (TTTj) . one has 

H=^-2H 2 , (13) 
o 



which can be used to express G from eq. (jl2|) as follows: 

G = AH 2 (R — 6H 2 ), (14) 

from which 

G = \HR 2 + 192H 5 - 32RH 3 + AH 2 R. (15) 
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Using 

^f(R,G) = f R R + f G G, (16) 
and ea . (| 1 5 p . R can be extracted from eq. (fTU|) as follows: 

(R - QH 2 )F R + GFq -F- 288H 2 (R/6 - 2H 2 ) 2 (F RG + AH 2 F GG ) 



GH(F RR + 8H 2 F RG + 16ff 4 F GG ) 

(17) 

H=^-2H 2 , (18) 



R 



where the second equation is the same as eq. (fT5|) . The set of above equations are 
the autonomous equations of F(R, G) gravity. The phase space of this problem 
is the two-dimensional (R — H) space. In the right-hand side of eq. (Tr7f . the 
expression (fTi]) must be used for the Gauss-Bonnet invariant G. Therefore the 
above equations are in the form 

H = fi(R, H), 

. (19) 
R = f 2 (R,H). 

The critical points are found by setting the eqs. ([T7|) and (|18p equal to zero. The 
result is 

^RF R + GF G - F = 0, (20) 
R=12H 2 . (21) 

Equation (TH|) also results in 

E>2 

G = 24H 4 = — (22) 
6 v ; 

at critical points. In obtaining eq. (|20|) . it is assumed that the denominator of 
eq. (|17|) has finite value at critical points. We will return to this assumption in 
section 5. Note that in eq. ([2H)l . G must be replaced by eq. (|2"2")l . In the case 
of f(R) gravity i.e. F(R,G) = /(i?)/2« 2 , the eqs.^)-^ axe reduced to the 
corresponding relations in ref. [15]. 

Since the effective equation of state parameter is defined through 

«-=jr-i-§^. (23) 

at critical points where H = 0, one has 

Wcff — » lu c = -1 (24) 

which is a characteristic of de-Sitter space-time. 

To study the stability of each critical point, one must evaluate the eigenvalues 
of matrix 



dh/dH dh/OR 



R=Rc ,H=H c 



R c and H c denote the values of R and H at the considered critical point. The 
critical point is a stable attractor only when the real parts of all the eigenvalues 
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of matrix M are negative. For negative real eigenvalues, the stable critical point 
is called a node, and for complex eigenvalues with negative real parts, the stable 
attractor is called a spiral. 

For autonomous equations (fT7|) and (fTH|) . the matrix M becomes 

M=M" , (26) 



-2F R /A H 



R=R C ,H=H C 



where 

A = F RR + 8H 2 F rg + 16H 4 F GG . (27) 
Therefore the eigenvalues are 



A L2 = i -3H±J(3Hf-J^-4H^j (2S) 

L V ^ ' J R=R C ,H=H C 

It is clear that the real part of eigenvalues Ai and A2 are negative, if and only if 



' 3A 



> 0. (29) 

R=R C ,H=H C 



This is the condition of stability of the attractors of F(R, G) gravity. The 
attractors are node if (3H C ) 2 > Aij and are spiral if (2>H C ) 2 < A-q . 
An interesting observation is that for i?-independent Lagrangian 

F(R,G)=F(G), (30) 

T) = —AH 2 < and therefore Ai > 0. So all the critical points of these models 
are unstable. 



3 Some examples of stable attractors 

In this section we will discuss two classes of F(R, G) models which lead to stable 
attractors. 

3.1 F(R, G) = F(RG) models 

For these models, the equation (|2T)|) results in 

^xF'(x) - F(x) - 0, (31) 

where x = RG and F'(x) = -^F(x). This relation determines the critical values 
x = x c . The stability condition (f2"9")l leads to 



F,( - X) - 1 > (32) 



9xF"{x) +4F'(x) 



at x — x c - For F(x) = x n cases, it can be shown that eq. ((3"Tj) does not have a 
nontrivial solution, except for a very special case which we will discuss it later. 
For F(x) — x n — c, where c is a positive constant, eqs.([3"Tj) and |3^|) result in: 



1 - 3n/2 



l/n 

(33) 



6 



and 



10 12 

18 18' 



(34) 



respectively. As an explicit example, we consider n = 11/18 and c = 1: 

F(RG) = (i?G) 11/18 - 1. (35) 

Eq. ([33| then results in 



So 



x c = R c G c = ^ = (12) 18 /n. 
o 



(R C ,H C ) = (7.047,0.766). 



(36) 



(37) 



Numerical calculation of eqs.(TT0 |) -(fT2 |) results in Figs. (1)- (4) for phase space, 
R(t),H(t), and u>(t) behaviors, respectively. These figures show that the point 
([3T)l is a stable attractor of spiral type. 




Figure 1: The spiral paths in R- H plane of F(R, G) = (RG) 11 / 18 - 1 model. 



H 




Figure 2: The plot of H(t) of F(R, G) = (i?G) n / 18 - 1 model. 
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Figure 3: The plot of R(t) of F(R, G) = (RG) 11 / 18 - 1 model. 




Figure 4: The plot of u(t) of F(R, G) = (RG) 11 / 18 - 1 model. 



3.2 F(R, G) = R + Gf{R) models 

For this functional form of F(R,G), the critical point equation ((20|) . using 
eq . (|22| . results in 

R 2 f(R) = 6. (38) 
This equation specifies R c . The stability condition (|29p reduces to 



i? 3 /"(i?)/12 + 2 



> 1, 



(39) 



which must be calculated at R = R c , For example for f(R) = mR n functions, 
where m and n are some constants, eqs. (13"8l and (15^1) result in 



R,= 



6 

mn 



and 



3 < n < -1, 



(40) 
(41) 



respectively. Since n is a negative number, m must be chosen negative so that 
R c in (|40p becomes a real positive number. As a specific example, we consider 
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m = — 1 and n — —2, or 

F{R,G)=R-G/R 2 . (42) 

R c and H c then become 

(R c ,H c ) = (l±). (43) 

Numerical results for Lagrangian f4"2"|) are given by Figs. (5) and (6), which are 
the paths in R — H plane and u)(t), respectively. Fig. (5) verifies that the critical 
point ([43]) is a stable attractor. 




0.1 - 

Figure 5: The spiral paths of F(R, G) = R - G/R 2 model in R - H plane. 




Figure 6: The plot of w{t) for F(R, G) = R - G/R 2 model. 



4 Critical points with zero eigenvalue 

As was mentioned previously, an attractor is stable if, and only if, the real 
parts of all eigenvalues are negative. Now what we can say if one, or more, of 
the eigenvalues are zero. The appearance of zero eigenvalues may be rooted in 
nonindependency of the chosen dynamical variables. For instance if we consider 
the phase space of F(R, G) models by three-dimensional (R, H, G) space, one 



9 



can show that besides the two eigenvalues Ai and A2 of eq. (P5|) . we have a third 
eigenvalue A3 = 0. This indicates that we can reduce the dimensionality of our 
phase space. 

For the cases where the dimension of the phase space can not be reduced, 
the appearance of zero eigenvalue means that the linear approximation, which 
leads to the matrix (|25[) , is not adequate and we must consider the higher order 
approximations to determine the behavior of the critical points. See for example 
[17]. 

In our F(R, G) models, for the cases where 



' 3A 



- 0, (44) 

R=R a ,H=H c 



one has 

Ai = , A 2 = -3H, (45) 

and the higher order approximations must be used to determine whether the 
considered critical point is stable or not. This is done by a standard method 
which has been discussed, for example, in [17] . 

At first, the eigenvectors of matrix M for the eigenvalues Ai = and A 2 = 
— 3H must be calculated, which result in 

^=(Iatt) > L V (46) 



x 24H C J ' rz y 6H L 

respectively. Using the transformation matrix T = (ipt,ip2), the new phase- 
space basis (U, V), i.e. the normal basis, can be found from (R, H) as follows 

which results in 



18 f c (48) 

v =m: i2AHcH - R) - 

To translate the critical point from (R c , H c ) to the origin of the phase space, 
we introduce TZ = R — R c and H = H — H c . 1Z and H, up to second order, then 
become 

K = H C K - 2AH 2 C H + DK 2 + BH 2 + CRH, 
■ 1 2 ( 49 ) 

n = -n- ah c h - m 2 . 

6 

In above equations, the eas. lfTT)) and P^|) have been used and the coefficients 
D, B and C are defined by 

1 f d 2 f 2 (R,H) \ 1 f d 2 f 2 (R,H) \ f d 2 f 2 (R,H) \ 

D= 2{ OR 2 ) ' 2 V dm ) ' C ={ dRdH J' (50) 

f2{B, G) is one introduced in eq. ([T^|) and all derivatives are calculated at critical 
point (R C ,H C ). Note that the linear terms in the right-hand-side of eq. 



10 



result in the matrix elements of M in eq. (126[) for the case where 77 = 0, or 
-2F R /A = -24H 2 . 

For the new phase space, introducing U = U — U c and V = V — V c , the 
eqs.gHJ) and (03) then result in 

U= (f+ j^f) + v f + 2DH c (iU + V) 2 + j(U + V)(4W + V) , (51) 

1> = -U + ^) (^ + V) 2 -2^i/ c (4W + V) 2 -^(W + V)(4W + V) . (52) 

Now taking V = /i(Z/) = all 2 + bU 3 + • • • , the coefficients a and 6 can be found 
by using the chain rule 

V = h\U)U . (53) 



By this way, the problem effectively becomes one-dimensional. Using eas.(|5ip 
and (|52[) . the coefficients of W 2 -terms of eq. (j53"|) result in the parameter a as 
follows: 

a =--4- ( ~+32H c D + -^— + ^c) . (54) 



3H C \3 18H C 3 

The coefficient b can be also found by the U 3 terms. Using the expansion 
V = aU 2 + eq.([5"T1) leads to 



u =(l +32H < D+ ik + l c ) u2 + {l + 16HcD + ir^l c ) aU3 + --- 

= aU 2 + f3U 3 + ■■■ . 

(55) 



So for the cases with zero eigenvalues, the higher order terms, through eq. (|55|) . 
must be considered in studying the stability behavior of attractors. Note that 
the absence of the linear terms in eq.(l5"5"|) reflects the fact that Ai is zero and 
we must focus on the next-leading terms. The attractor is then stable if a < 0. 
For the cases where a = 0, we must look at the sign j3. (3 < leads to stable 
attractors. If again (3 becomes zero, we must go to higher orders. 
As an explicit example, we consider 

F(R,G) = R+±R 3 -± . (56) 

The critical points are obtained by solving eqs. ([20|) and (|2"Tj) . which result in 

(R C ,H C )= (1,4=) ■ (57) 
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Note that we have not considered the unphysical solution R c = — 2. It can be 
easily seen that 77 = (using ((2l)|) ), therefore Ai = and A2 = —2>H. Despite a 
negative eigenvalue, the numerical calculations show that the critical point (|57p 
is not a stable attractor (see Fig. 7) 

This can be justified by calculating the parameter a and [3 of eq. (|55j) . which 
leads to 

a = lw-™ & + .... (58) 

3 9V3 

Since a > 0, the critical point is not stable, in accordance with Fig. 7. 
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R 




10 20 30 40 50 

Figure 7: The plot of R(t) of F(R, G) = R + ±i? 3 - §. The initial values are 
(Ro,H ) = (1.1,0.3). It is clear that system does not approach R c = 1. 



5 Singular points of F(R, G) 

As pointed out after eq. . in deriving the critical point equation ([20]) . it has 
been assumed that the denominator of eq. (fT"7|) is finite at the critical values, so 
R = leads us to set the numerator of eq. (fTT)) equal to zero. But, as we will 
show, for the cases where the Lagrangian F(R, G) has some singularities, this 
assumption, i.e. the finiteness of the denominator of ea. (|17|) . is not right and 
we must carefully reinvestigate our results. 
Take the function F(R, G) as follows 

F < R - G >=§t§ < 59 ' 

where Q(R, G) has a root of order n at R — a, i.e. 

Q(R,G) = (R-a) n g(R,G). (60) 

It can be easily seen that in this case, the denominator of eq. (fTT|) diverges at 
R — > a, with the power greater than numerator. Substituting cq. (|59l) into 
cq.fHU), results in 

■ = (6H 2 - R)PQQ R + 8H 2 {R - 12H 2 f(2PQ R Q G - PQQrg - QPgQr) + ■ ■ 

6H(2Q 2 R - QQ RR )P + ■■■ 

^ (61) 

where the dots denote the higher order terms of R — a. Power counting of 
eq. ([6"Tj) shows that 

■ Q({R-af^) 

R ~ 0((R-a)*"-*) ' (62) 
which results an extra solution for R = as follows 

(Rc,H c )= ■ (63) 

H c is found from eq. (I^T|) . 
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To study the stability of this critical point, we need to calculate the eigen- 
values of the matrix M (in eq . (j2l)|) ) . A lengthy calculation shows that 



SR 



-H 



Qr + QQrr QQr(3QrQrr — QQrrr) 



2Q 2 r -QQrr (2Q 2 r -QQrr) 2 
where using the expression ([6"U|) for Q(R, G), leads to 



SR 



SR= - 



H, 



at R 



n + l 

a. So the matrix M becomes 



with eigenvalues 



M 



Ai 



-4# c 




SR + --- , 



-1/6 
-H c /(n+l) 



A, 



n + l 



(64) 



(65) 



(66) 



(67) 



where both of them are real negative numbers. So we lead to an important gen- 
eral consequence: Any singularity of the function F(R, G) is an stable attractor 
solution. 

The same is true for the cases where Q(R, G) has a root of order n at G — j3: 



Q(R,G) = (G- P) n g{R,G) . 
The same procedure results in a critical point at G c — f3, or 

l/4\ 



{R C ,H C ) 



(68) 



(69) 



in which eq. ([^f has been used. The matrix M becomes the same as eq. (l6l)l) . 
which proves that this critical point is a stable attractor. 
So generally for 



Q(R, G) = (R- a) n (G - (3) m g(R, G) 



(70) 



the model has the stable attractor points (|6"3"|) and (|6"9"|) . 

As an example, we consider the model discussed in section 3-2, that is 
F(R,G) = R + Gf(R). The regular critical points (non-singular type) can 
be found by solving the relation (l38|) . and the stability condition is eq. ([39|) . 
Now consider the explicit example 



F(R,G) = R 



where m is a constant. Eq. 



3(R 2 - l) 2 



leads to 

mR 3 = 

and the inequality (|39[) results in 

2(fl c 2 -l) 
R2-5 



mG 
R^~T ' 



R c = R c {m) , 



> 1 



(71) 



(72) 



(73) 
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For i?;? — 5 > 0, eg . ([75]) leads to fZj? > —3, which is always true, and for 
Rl — 5 < 0, it results in R? c < —3, which is never true. So the condition ([73]) 
holds if 

R 2 C > 5 . (74) 
Now if we chose m = — 10, the eq. (|72| gives two following real solutions: 

R lc = 0.534 , R 2c = 3.837 . (75) 

It is clear that R\ c does not satisfy (l?4"|) . while i?2c does. Explicit calculation of 
rj (in eq.(Jll|)) shows that ti\r=r 1c < and r?|ij=ii 2c > 0. So we expect that the 
stable critical point of 



F{R,G) =R 



IPG 
R?-l ' 



is 



(i? c ,# c ) = (3.837,0.565) 
Numerical calculation verifies this. See Fig. 8 



(76) 
(77) 



4.4 




4.2 




4.0 




3.8 




3.6 




3.4/ 




3.0 





Figure 8: The spiral paths leading to attractor (jTTJ) of the Lagrangian 
F(R,G) = R-10G/(R 2 - 1). 

Until now, we find the regular attractor of (f76|) . But it is clear that the 
F(R,G) in eq. ([76|) is singular at iZ = 1. So we expect another stable attractor 
at the point 

(Rc,H c )= ■ (78) 



This new attractor is also verified by the numerical method. See Fig. 9 



6 The critical curves 

There are other interesting cases in which the critical points are replaced by 
critical curves. In this case, each of the infinite points on this critical curve are 
in fact a critical point, and besides, as we will show, they are stable attractors. 
This situation occurs when the criticality condition holds for any R and H 
values. 
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0.70 - 

Figure 9: The paths of F(R, G) = R - WG/(R 2 - 1) leading to attractor 



Before introducing some special examples, we first prove a general statement: 
If a F(R, G) function satisfies (20]) and R and G satisfy eqs.j2T|) and ([22]). 
respectively, then r] = Fr/(3A) — AH 2 is equal to zero. 

proof: Since condition (|2"0"f holds for any R and H, it can be differentiated with 
the result 

{\rFrr - ^F R + GF RG )dR + (^RF RG + GF GG )dG = 0. (79) 

But from eq. ([T4"|) we have 

dG = 8H(R~12H 2 )dH + 4H 2 dR, (80) 

so 

[\rFrr - \f r + GF rg + 4H 2 (IrF rg + GF GG )]dR 

x (81) 
+ 8(-RF RG + GF GG )H(R - 12H 2 )dH = 0. 

Using eq. (|2"Tj) , the coefficient of dH becomes zero. The coefficient of dR, which 
now must be set to zero, specifies Fr as follows 

F R - R{Frr + 8H 2 F RG + 16H 4 F GG ) = RA, (82) 

in which the function A, introduced in (127p . has been used. Therefore r\ becomes 

Fb 1 
V = ^j-4H 2 = -R-4H 2 = 0, (83) 

where R = 12H 2 has been used. This completes our proof.B 

Now we note that if F(R, G) satisfies (|20p . this equation does not impose any 
extra constraint on F(R, G) and therefore does not specify any critical values 
for R and H. So the only remaining relation in the R — H phase-space plane is 
the second equation (f2Tj) . which defines a critical curve. The eigenvalues of this 
critical curve, as a result of the above statement which leads to eq. (|83p . are : 

Ai = , A 2 = -3H. (84) 
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But it can be shown that in the case of the emergence of critical curve, the 
stability of any particular point on this curve can be determined by the nonzero 
eigenvalues [15]. Since in our case, A2 = — 3H < 0, therefore any points on 
the critical curve R = 12H 2 of F(R,G) models is a stable attractor. This is a 
general result. 

Now let us consider some explicit examples. 

Example 1: Let us first consider the class of models introduced in section 
3-1. For F(R, G) = F(RG), it is obtained that eq.^Q results in 

^xF'(x) - F{x) = 0, (85) 

where x — RG. If we demand that the eq. ([85| satisfies for all x, then it can 
be viewed as a differential equation with solution F(x) = a; 2 / 3 . So the F(R, G) 
model 

F(R,G) = (RG) 2/3 (86) 

has a critical curve R = 12H 2 . All the points on this curve are stable attractors. 
Figs. (10) and (11) show that both H lc = 0.977 and H 2c = 1-99 points (as two 
arbitrary points), with Ri c = 11.454 and R 2c = 47.52, respectively, are stable 
attractors of this model. 

H 

0.9785 f 



0.9780 
0.9775 



2 ^— 4 6 8 tO 



Figure 10: The plot of H(t) for F(R,G) = [RG) 2 / 3 model. The point 
(H c , R c ) = (0.977, 11.454) is a stable attractor. 

Example 2: consider the following F(R, G) model: 

F(R, G) = R n f(G k R m ) (87) 
with arbitrary constants n, k and m. Substituting ()87f) into eq. (|20p . results in: 

(| - l)R n f{x) + (~ + k)R n+m G k f'(x) = (88) 

where x = G k R m . If one demands the above equation satisfies for all R and Gs 
and for any arbitrary function f(x), then the constants n, k and m must satisfy 

n = 2 , m = -2k. (89) 

So F(R, G) = R 2 f(G k /R 2k ) = R 2 g(G/R 2 ) satisfies (J20j) and the curve R = 
12H 2 is its critical curve. It is interesting to note that the case considered in 
example 1, i.e. the eq.®, is in fact of this type: (RG) 2 / 3 = R 2 (G / R 2 ) 2 ' 3 . 
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Figure 11: The plot of H(t) for F(R, G) = (RG) 2 / 3 model. The point 
(H C ,R C ) = (1.99,47.52) is a stable attractor. 

Example 3: Consider the following model 

F(R, G) — uG + f(R) , (90) 

then eq. (j2"U|) results in: 

Rf = 2/ (91) 

which its solution, as a differential equation, is f(R) — (3R 2 . So all models of 
the type 

F(R, G) = aG + [3R 2 , (92) 

have R — 12H 2 as their critical curve, with infinite number of stable attractors. 

This example also shows that in the case of f(R) gravity theories, which 
is the a = case of eq. ([90|) . the only model which leads to the critical curve 
R = 12H 2 is the f(R) = R 2 model. 

Example 4: As the last example, consider the model 

F(R,G) = R + Gf(R). (93) 
Substituting (f9"3f into eq. ([2TJ]) , results in 

R 2 f'(R)=6^ /W = -|- (94) 

So F(R, G)=R- 6G/R has also the critical curve R = 12H 2 . 

The above mentioned procedure can be applied to some other functional 
forms, such as F(R, G) = f(G/R), with the result F = (G/i?) 2 , etc. 

7 Conclusion 

As a candidate of dark energy, we consider the generalized Gauss-Bonnet dark 
energy models looking for the situations where the late-time behavior of this 
modified gravity theory is the de-Sitter space-time. We describe the phase 
space of this theory by the two-dimensional R — H space. This dimensionality 
verifies by the fact if the three-dimensional R — H — G space has been chosen, 
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one of the eigenvalues of stability matrix is always zero, which indicates that 
the number of independent variables is two. 

The eigenvalues of stability matrix show that the critical points, i.e. the de- 
Sitter space-times, are, in general, the stable attractor if t] = Fr/(3A) — AH 2 > 
0, a fact that has been verified by several examples. The emergence of critical 
points with rj = 0, in which one of the eigenvalues is zero, Ai = 0, forces us 
to consider the higher order terms in normal basis in order to have a correct 
judgment about the stability of these kinds of critical points. 

We also find two classes of stable attractors: the singular points of the 
Lagrangian F(R, G) and the cases where the critical points are replaced by the 
critical curve R = 12H 2 (in R — H plane). In the latter case, all the points on 
this curve are stable attractors. 
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